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The NJL-jet model provides a framework for calculating fragmentation functions without the 
introduction of ad hoc parameters. We develop the NJL-jet model to investigate dihadron frag- 
mentation functions (DFFs) of the form D^ 1 ' h2 (zi, Z2). Here we studied DFFs for q — > {ir + ir~}, 
{tt + K~} and {K + K~} with q = u,d, s. The driving terms, which represent the probability of one of 
the hadrons being emitted in the first emission step of the quark-jet hadronization picture, dominate 
the solutions of the DFFs where either z\ or Z2 is large, and z\ (22) is the light-cone momentum 
fraction of the emitted hadron, h\ ihz). The higher order terms, which represent the probability 
of neither of the hadrons being emitted in the first emission step of the quark-jet, become more 
significant as z\ (22) is lowered. Finally, we present a sample result for QCD evolution of DFFs, 
that significantly modify the model solutions when evolved to typical hadronic scale of 4 GeV 2 . 
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I. INTRODUCTION 

Deep inelastic scattering (DIS) has proven to be an 
invaluable source of information about the structure of 
the nucleon pQ. Initially it provided critical information 
on the relative distribution of momentum between va- 
lence and sea quarks and the gluons. As the experimental 
capabilities have grown so have our ambitions and over 
the past decade semi-inclusive deep-inelastic scattering 
(SIDIS) has helped (along with Drell-Yan) to expand our 
knowledge of flavor structure and asymmetry [2HZ] ■ With 
several new experimental facilities with 100% duty factor 
under construction, SIDIS will become even more impor- 
tant. For example, we may finally be able to pin down the 
elusive s—s asymmetry [8llll|. Another area of great cur- 
rent excitement concerns the distribution of the spin of 
the proton [12Tl32j . There polarized SIDIS is potentially 
extremely valuable through the study of TMD's [551142) . 
which will complement work on GPD's [H)l |4"5H4"B] . 

For these studies to achieve their full potential it is 
vital that we develop the deepest understanding of the 
fragmentation functions [39], especially their flavour de- 
pendence and ultimately their dependence on spin and 
transverse momentum. Fragmentation functions are an 
important theoretical tool in the investigation of scat- 
tering reactions, for example in the separation of the 
flavor dependence of parton distribution functions. Ex- 
perimental extractions of fragmentation functions from 
deep-inelastic scattering data [50] [5TJ have increased the- 
oretical activity in this area. Yet the phenomenological 
extraction of even favored fragmentation functions suf- 
fers from significant uncertainty while the situation for 
the unfavored is worse. This in turn effects the system- 
atic errors associated with the extraction of the flavor de- 
pendence of parton distribution functions through SIDIS. 
These considerations have led us to develop and study the 
NJL-jet model [521455] , This model builds on the Field- 



Feynman Quark-jet model (FFQJM) [56], by using an 
effective chiral quark model to provide a framework in 
which calculations of both quark distribution and frag- 
mentation functions can be performed without introduc- 
ing ad hoc parameters. Pion fragmentation functions in 
the NJL-jet model were calculated in Ref. [55]. The NJL- 
jet model was then extended to include strange quark 
contributions and kaon fragmentation functions were ob- 
tained [53 . Further extensions of the model are the in- 
clusion of vector meson, nucleon and antinucleon frag- 
mentation channels [54 , and the inclusion of transverse 
momentum dependence [55 . 

Dihadron fragmentation functions (DFFs) represent 
the probability of producing two hadrons in the decay 
chain of a fragmenting quark. Some recent work in the 
area of DFFs include Refs. [57] [55]. In Ref. [57], param- 
eters for a spectator model are fitted to output from the 
PYTHIA event generator [59] tuned for HERMES [60] 
for dihadron fragmentation functions with a dependence 
on the sum of the light-cone momentum fractions of the 
two produced hadrons and their invariant mass squared. 
DFFs for large invariant mass are studied in Ref. [58] . 
The dihadron fragmentation functions evolution equa- 
tions are derived in Ref. |61) from factorization of the 
cross-section for the production of two hadrons in e + e~ 
annihilation in the MS factorization scheme. Evolu- 
tion equations for non-singlet quark DFFs are studied 
in Ref. [B2], while the singlet quark and gluon DFF evo- 
lution equations are studied in Ref. [53]. In Refs. [6"21 l6"3l 
the ratio of the dihadron and single hadron fragmenta- 
tion functions are examined, as this ratio is useful when 
considering experimental measurements. The choice of 
initial conditions is studied in Ref. |64j , primarily by con- 
sidering the two-body correlation function. 

Depending on the polarization of the fragmenting 
quark, special types of DFFs known as Interference frag- 
mentation functions (IFFs) can be constructed. The 
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chiral-odd IFFs can be related to transversity [qT5TI69| . 
Transversity is one area of current interest that requires 
knowledge of the fragmentation functions of quarks [27J . 
Out of the three leading-twist distribution functions that 
describe the quark structure of hadrons, it is the least 
well known, the other two being the unpolarized and hc- 
licity distributions. Recent work [T0HT2"] suggests that di- 
hadron fragmentation functions (DFFs) may be useful in 
extracting transversity distributions by considering the 
semi-inclusive deep-inelastic scattering (SIDIS) produc- 
tion of two hadrons with small invariant mass. Though 
transversity will not be the focus of this paper, it is pre- 
sented as one possible motivation for further investigation 
into DFFs. 

In this work we extend the latest version of the NJL- 
jct model, presented in Refs. [531 [M], to investigate di- 
hadron fragmentation functions. In Section[ll]we present 
a summary of fragmentation functions in the NJL-jet 
model, as set out in the aforementioned papers, with a 
focus on those parts that are relevant to understanding 
the dihadron fragmentation functions. Section III out- 



lines the extension of the NJL-jet model to be used in 
investigating DFFs, while results at the model scale for 
the DFFs are presented in Section IV In Section [V] we 
briefly discuss the QCD evolution equations for DFFs 
and present sample evolution results for our model. 



II. QUARK FRAGMENTATION FUNCTIONS 

This section provides a quick overview of the calcula- 
tion of the quark fragmentation functions in the NJL- 



jet model 52~{54| . focusing on the aspects important 
to obtaining dihadron fragmentation functions within 
the model. Here we employ the 5/7(3) Nambu-Jona- 
Lasinio (NJL) effective quark model [TSTfT?] using light- 
cone (LC) coordinates [53] . In the NJL model we include 
only the four-point quark interaction in the Lagrangian, 
with up, down, and strange quarks and no additional 
free parameters. We employ Lepage-Brodsky (LB) "in- 
variant mass" cut-off regularization for the loop integrals 
(see Refs. [53] for a detailed description as applied to the 
NJL-jet model). 
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FIG. 1. Cut diagrams for quark fragmentation function. Solid 
lines represent quarks and the double dashed lines a meson. 



The quark fragmentation function D^{z) is the prob- 
ability for a quark of type q to emit a hadron of type h 
carrying fraction z of its light-cone momentum h (here 
meson m — qQ). We denote the elementary quark frag- 
mentation function, corresponding to the situation where 
the detected hadron is the only emitted hadron, by (z). 
The corresponding cut diagram for the elementary quark 
fragmentation function is shown in Fig. [T] The elemen- 
tary fragmentation function depicted in Fig. [T] can be 
written as: 
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where C™ is the corresponding flavor factor and g m qQ is where 
the quark- meson coupling. The masses Mi, M 2 and m m 
are the masses of the fragmenting quark, the remnant 
quark and the produced hadron (here the hadron is a 
meson), respectively. 
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If a sharp cut-off in the transverse momentum 
is assumed, the integration in Eq. ([!]) can be evaluated 
analytically (Eq. S). 



A=({z~l)M 1 +M 2 ) 2 , 

B = z(z- l)Ml + zMl + (I - z)m 2 m 



(3) 
(4) 



The Lepage-Brodsky "invariant mass" cut-off regular- 
ization method (Refs. [55] and [75] describe this when 
applied to the NJL- Jet model) is employed to regularize 
the loop integrals. The loop integrals are regularized by 
setting a cut-off on the invariant mass, Mi 2 , such that 
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where Ai 2 is the maximum invariant mass. Here the 3- 
momentum cut-off, denoted by A 3 , is fixed by reproduc- 
ing the value of the experimentally measured pion decay 
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FIG. 2. Quark Cascade 



constant. In Lepage-Brodsky regularization, is given 
by: 



Pi 
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The value of the 3-momentum cut-off used in this 
work, A 3 = 0.67 GeV, was obtained in Ref. [75] by 
choosing the constituent light quark mass M = 0.3 GeV 
and using pion decay = 0.093. The correspond- 
ing constituent strange quark mass, M s — 0.537 GeV, 
was obtained by reproducing the pion and kaon masses, 
= 0.14 GeV and m K = 0.495 GeV. The calculated 
values of of the quark-meson couplings are g^qQ = 3.15 
and gKqQ = 3.3876. 

The elementary quark fragmentation function, inte- 
grated over the light-cone momentum fraction z, rep- 
resents the total probability of a quark splitting into a 
hadron of a given type plus another quark. In construc- 
tion of the NJL-jet model we are interested in procesess 
where a hadron is produced at each step. Thus we con- 
struct the renormalized elementary fragmentation func- 
tions, d 1 q n (z), such that the total probability of emitting 
a hadron (summed over all possible hadron types m' that 
the quark q can emit in the elementary splitting process) 
is one: 
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In a quark-jet-model, the total fragmentation function, 
D q n (z), is described by successive elementary splittings 
of a quark into hadrons. This process is shown diagram- 
matically in Fig. [2j The initial quark q fragments into 
a meson, m = qQ, with light-cone momentum fraction 
z of the initial quarks momentum, and a quark, Q, with 
light-cone momentum fraction 1 — z. The emitted quark 
Q fragments as well, and the process repeats, forming a 
cascade of hadrons. It is important to note that within 
the model the emitted hadrons do not interact with the 
other hadrons produced in the quark jet. An integral 
equation for the quark cascade process shown in Fig. [2] 
was derived in the quark-jet model of Ref. |56j . The in- 
tegral equation for the total fragmentation function is: 
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dy 
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where df(z) = d™(l - z)\ m=qQ . 

The probabilistic interpretation of the Eq. ^ can be 
clarified by multiplying both sides by a factor of dz. The 
term on the left hand side is the probability for the quark 
q to emit meson m with light-cone momentum fraction 
z. On the right hand side, the first term is the driving 
function, which represents the probability of creating a 
meson m carrying momentum fraction z to z + dz from 
the first emission step and the second term represents the 
probability of creating the meson, m, further in the quark 
decay chain. The above equation is solved by uniformly 
discretizing z and y in the interval [0, 1] and approximat- 
ing the integrals as sums over these discrete values of z 
and y. Then D™(z) and d™{z) can be expressed as vec- 
tors D™ and / of values at the discretization points of 
z and the integrand of the second term, without D™{y), 
can be written as a matrix g over the values of the dis- 
cretization points of z and y: 



D q n = f + g-D™ 
(!-<?)■ D™ = f 



D™ = (I-g)- 1 f, 



(9) 



where / is the unit matrix. 

Here it is important to use an appropriate number of 
discretization points to avoid large numerical errors when 
solving for D q 1l {z). The number of points used was in- 
creased until there was sufficient convergence of the solu- 
tions of the fragmentation functions. The resulting solu- 
tions for the fragmentation functions of u, d and s quarks 
to 7T + and K + are presented in Ref. [54] . The fragmenta- 
tion functions of u, d and s quarks to n + are also shown 
here in Fig. |3] as well as the driving function of the u 
to 7r + , d^ , which is labelled as Df : u (the notation Df : q 
will be adopted in Section IV A as well). 
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FIG. 3. Fragmentation Functions of u (black solid line), 
d (blue dashed line) and s (red dot-dashed line) quarks to tv + 
as a function of the light-cone momentum fraction z. Driving 
function of the u (purple dotted line) to 7r + is labelled as Df:u. 
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III. DIHADRON FRAGMENTATION 
FUNCTIONS 

We now consider a semi-inclusive process in which 
two hadrons are detected in the final state. This re- 
quires a new fragmentation function, known as the di- 
hadron fragmentation function, that describes the prob- 
ability of this process. Dihadron fragmentation func- 
tions (DFFs) may be useful in the extraction of transver- 
sity distributions [70], which are the least well known 
of the three leading-twist distribution functions that de- 
scribe the quark structure of hadrons, the other two be- 
ing the unpolarized and helicity distribution functions. 
We now extend the NJL-jet model to describe the DFFs. 
Dihadron fragmentation functions, D^ 1 ' /t2 (z 1 , z 2 ) 1 corre- 
spond to the probability of a quark q producing two 
hadrons, hi and h 2 , that carry its light-cone momentum 
fractions Z\ and z 2 , respectively. An illustration of how 
a quark cascade can produce two observed hadrons, hi 



and h 2 , in the NJL-jet model is shown in Fig. [2j 

The integral equation for the dihadron fragmentation 
function D^ 1,h2 (zi, z 2 ) has been constructed by Field and 
Feynman in the quark-jet model (Eqs (2.43a)-(2.43d) of 
Ref. [56]), which is shown in Eq. (10). Here d^(z) and 



dqiv) are the elementary splitting functions of the quark 
q to the corresponding hadron h and quark Q. On the 
left hand side of Eq. ( 10 ) is the term representing the 
probability for the quark q to emit hadrons hi and h 2 
with light-cone momentum fractions Z\ and z 2 , respec- 
tively. The first term on the right hand side of Eq. ( 10 ) 



corresponds to the probability of producing hadron hi 
from the quark q at the first step in the cascade, followed 
by hadron h 2 produced either directly afterwards or fur- 
ther down in the quark decay chain, while the second 
term is similar to the first one, except for hi •<-> h 2 • The 
third term corresponds to the probability of having both 
hi and h 2 produced after the first hadron emission. 



D^ h Hzi,z 2 )=d h ^(zi) 



rth 2 ( _Z2_ 

^qi \ l-zi 
l-zi 



1-22 



Q J Zl+22 ' 



V V 



(10) 



q^hi+qi: q^h 2 + q 2 . 
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In the integral term we perform a change of integra- 
tion variables to £i = Zi/rj and £2 = z 2 /rj, so that the 



arguments of D'q'* 12 (£1, £2) will correspond to £1 and £2 
at grid point values when uniformly discretized: 
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Then the equation for the dihadron fragmentation functions takes the following form: 

D h 2( z 2 ) D hU 21 \ 

D^( Zl ,z 2 ) = {zi ) + ) -^T^ 2 
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To solve the above equation for the dihadron fragmen- the integrals as sums over the discretized values of these 
tation function -D? 1 ' 2 {zi,z 2 ), we discretize Zi, z 2 , £1, variables. The fragmentation functions are written in 
and £2 uniformly in the interval [0, 1] and approximate matrix form, where the elements of the matrices are their 
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values at the corresponding uniformly discrete values of 
the arguments. We used Mathematica to solve for both 
the single hadron and dihadron fragmentation functions. 
The number of discretization points used for the single 
hadron fragmentation functions was 500, while the num- 
ber of the discretization points afforded for the dihadron 
fragmentation functions was 200. These values for the 
numbers of discretization points produced convergence 
of the solutions within typically 5%, while allowing for 
a reasonable computational time and computer memory 
size required by the problem. Several techniques were 
used to lower the memory use of the program, including 
the use of the sparse arrays in Mathematica. To calcu- 
late the third term of Eq. ( 13 1, the integrals over £i were 



converted to a sum over its uniformly discrete values. 
The delta function was used to eliminate the integration 
over £2- The values of £2 that are selected by the delta 
function may not match any of its uniformly discretized 
values. To account for this, the values of the DFFs at 
the selected £2 were obtained using linear interpolation 
from neighboring discrete values. 



IV. RESULTS 

In this section we investigate various features of the 
DFFs obtained as solutions of Eq. (fl3|) . Section |IV A 



investigates the contribution of the integral term versus 
that of the driving term to the solution of DFF. The 
impact of the inclusion of the strang e quark o n * 



IV B 



In Section 



IV C 



we 



and * is studied in Section 

consider D^ +K (21,22) an d fix either Z\ or z%, to study 
the dependence of this DFF on each of the variables. 



A. Contribution of the Integral Term 



In Eq. ( 13 ), the sum of the first two terms is considered 



to be the driving function of the dihadron fragmentation 
function and they describe the probability of emitting one 
of the detected hadrons in the first emission step. The 
last term in Eq. ( 13 ) corresponds to the the probability 



of emitting both detected hadrons after the emission of 
a hadron in the first step. We now consider the contri- 
bution of this last term to the solution for DFF by com- 
paring them with the corresponding driving functions for 
three combinations of observed pions and kaons: 7r + 7r~ , 
tt+K- and K+K-. 



The DFFs for the produced hadrons ir^ 



tt + K and 



K + K as functions of 2:2 with z\ = 0.5 are shown in 
Fig. |U The plots in the fi gures show that the favored 
DFFs, where the initial quark can produce either of the 
detected hadrons from the initial quark, are almost equal 
to the driving function, with the integral term giving only 
a very small contribution. The unfavored DFFs, where 
neither of the hadrons can be directly produced by the 
initial quark, are generated entirely by the integral term. 



In Figs. [4] and [5] the solution of the DFF for the up 
quark is shown by the orange circle points and the driving 
function is shown as a solid gray line. The green diamond 
points and dotted black line show the DFF and the driv- 
ing function for the down quark, respectively. The DFF 
and the driving function for the strange quark are shown 
by the blue square points and the dot-dashed red line, 
respectively. In each of the figures, the number in the 
brackets in the corresponding legend indicate the scaling 
factor used in depicting the curve on the plot. This no- 
tation for the scaling factor is also used in Sections |IVB| 
and UVCl 

In Fig. [5j the results for the DFFs and driving func- 
tions for the 7r + 7r~, ir + K~ and K + K~ are presented for 
the fixed value of z\ = 0.1. Here the integral term con- 
tribution to the up quark DFFs become visible as the 
value of z\ is lowered because the driving functions con- 
tribution to the DFF becomes less significant. The driv- 
ing functions contribution for the down (— > 7r + 7r~) and 
strange quark (— s- n + K ~ and K + K~) DFFs are still very 
dominant, so there is no noticeable contribution from the 
integral term here. It is worth noting that the integral 
term contributions to both favored and unfavored DFFs 
are of the same magnitude, but the contributions to the 
favored DFFs are only noticeable when the driving func- 
tion is not dominating the solution. 



B. Impact of Including the Strange Quark on the 

7"")7T It 
U 1 



We now study the impact of the inclusion of the strange 



and Dl * . In Eq. the integral 



quark on 

term contains a sum over Q that runs over the flavors 
of the quarks considered in the model, thus the inclusion 
of the strange quark couples the DFFs for the u and d 
quarks to those of the s quark. Also, the inclusion of 
the strange quark affects the single hadron fragmenta- 
tions of the driving terms, as in their respective integral 
equations there is a sum over Q as well (Eq. (|8|). This 
potentially can have a large effect, as in Section [IV A| it 
was shown that the driving functions give the most of the 
contribution to the favored DFFs. 

The solution of * for z\ fixed at 0.1 and 0.5 are 
shown in Fig. [6] Here, the dashed blue line and the red 
dot-dashed line represent the results for up quark DFF 
with and without the strange quark (denoted by (s) and 
(no s) in the captions), respectively. Similarly, the purple 
dotted line and the black solid line represents the results 
for the down quark DFF with and without the strange 
quark. 

The shapes of the dihadron fragmentation functions re- 
main the same for both the u and d quark DFFs, with the 
down quark DFF being the larger in magnitude compared 
to the up quark DFF for low z± and vice versa when z\ 
is increased. The main change is the considerable reduc- 
tion in the magnitude of DFFs when the strange quark is 
included, caused by the availability of the kaon emission 
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FIG. 4. Dihadron fra gme ntation functions for Z\ 



= 7T + ,h 2 

h 2 = K- 



7T , 

The 



0.5 for[(a)1 

le up, 



(b) hi = n + ,h 2 = K~ and 
i 2 = i\ . ± ne DFFs and driving functions of t 
down and strange quarks are shown by the orange circle points 
and solid gray line, green diamond points and dotted black 
line, and the blue square points and the dot-dashed red line, 
respectively. Driving functions for fragmenting quark q are 
also labelled in legend as Df:q. The number in the brackets in 
the legend indicates the scaling factor used in depicting the 
curve. 



FIG. 5. Dihadron fra gme ntation functions for Z\ = 0.1 f or[(a) 



hi 
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= 7T + ,h 2 

M = K- 



(c) 



hi 
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(b) hi = n + ,h 2 = K~ and 
The DFFs and driving functions of t 
down and strange quarks are shown by the orange circle points 
and solid gray line, green diamond points and dotted black 
line, and the blue square points and and the dot-dashed red 
line, respectively. Driving functions for fragmenting quark q 
are also labelled in legend as Df:q. The number in the brackets 
in the legend indicates the scaling factor used in depicting the 
curve. 
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channels. Thus the inclusion of the strange quark in our 
model proves to be very important in describing the light 
quark DFFs. 
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FIG. 6. Comparison of strange quark con tribution of 7r + 7r~ 
dihadron fragmentation functions for (a) z\ = 0.1 and (b) 
z\ — 0.5. The dashed blue line and the red dot-dashed line 
represent the results for up quark DFF with and without the 
strange quark (denoted by (s) and (no s) in the captions), 
respectively. Similarly, the purple dotted line and the black 
solid line represents the results for the down quark DFF with 
and without the strange quark. The number in the brackets 
in the legend indicate the scaling factor used in depicting the 
curve. 



C. Study of D\ K 

In this section, we examine the plots of D™ +K , where 
either z\ or z 2 is fixed. These particular DFFs were cho- 
sen since q — > tt + K~ is a favored fragmentation chan- 
nel to one of the hadrons both for a light and a strange 
quark q. This produces more interesting results to ex- 
amine than if we had chosen * , as those DFFs are 
symmetric in q = u and q = d, thus the DFF for q = u 
at fixed z\ is same as that for the DFF for q — d at 



fixed 22, etc. The results for fixed values of z\ and z 2 
are shown on the plots in Figs. [7] and [HJ The up quark 
DFFs are represented by dotted red lines, while the down 
quark DFFs are represented by dashed blue lines and the 
strange quark DFFs are represented by solid black lines. 
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FIG. 7. 7T K dihadron fragmentation functions for |(a)| zi = 
0.1 and (b) z 2 =0.1. The up, down and strange quark DFFs 



are represented by dotted red lines, dashed blue lines and solid 
black lines, respectively. The number in the brackets in the 
legend indicate the scaling factor used in depicting the curve. 



We first examine the DFFs for z\ = 0.1 (Fig. 7(a) ) and 
z 2 = 0.1 (Fig. |7(b)[ ). Since the hadron corresponding to 
the fixed light-cone momentum fraction only has a small 
amount of the fragmenting quarks momentum, most of 
the momentum is attributed to the favored fragmenta- 
tion channel of the other hadron. For Figs. |7(a) and |7(b){ 
this corresponds to the strange and up quarks fragmenta- 
tions to the K~ and 7r + , respectively. The down quark is 
unfavored for both hadrons and thus receives very little 
contribution to its DFF in both plots. 

After increasing the fixed value of Z\ (Fig. 8(a) ) and 
z 2 (Fig. |8(b)| to 0.5, the strange quark DFFs are the 
largest for both. The strange quark is the only initial 
quark that can produce both hadrons in the first two 
steps of the cascade, whereas, both the up and down 
quarks require multiple decays to produce both hadrons. 
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FIG. 8. ty + K dihadron fragmentation functions for |(a)| Z\ = 
0.5 and (b) 22 = 0.5. The up, down and strange quark DFFs 



are represented by dotted red lines, dashed blue lines and solid 
black lines, respectively. The number in the brackets in the 
legend indicate the scaling factor used in depicting the curve. 



The up quarks solutions are the second largest since it can 
produce the n + hadron in the first emission step, while 
the down quarks solutions are low for both plots, as it 
can't produce either of the hadrons in the first emission 
step. 



V. EVOLUTION OF THE DFFS 

The results for the dihadron fragmentation functions in 
the NJL-jet model presented in Section |Tv| are all at the 
model scale of 0.2 GeV 2 . The model scale was obtained 
in Ref. [S3] such that after NLO evolution the model 
calculated u quark valence distribution function in the 
7r + matched those experimentally measured in Refs. [75] 
and [80) . To compare our results to experiment or re- 
sults from other models, we need to evolve the dihadron 
fragmentation functions to an appropriate momentum 
scale. In Ref. [61] . the dihadron fragmentation func- 
tions evolution equations are derived from factorization 
of the cross-section for the production of two hadrons in 
e + e~ annihilation in the MS factorization scheme. Using 
Jet-Calculus, Ref. [5TJ deduces the evolution equations 
for DFFs with an explicit dependence on the invariant 
mass of the hadron pairs, Mh- The DFFs with a depen- 
dence on the invariant mass are addressed as extended 
dihadron fragmentation functions (extDiFF). The extD- 
iFFs are important as they will relate to experimental 
results that include the dependence on invariant mass 
spectra. 

The leading order (LO) evolution equation for DFF 
from Ref. [ST] is presented in Eq. (14). 



d 



dlnQ z 2n J,.i_,„ it J \u u J 



g s (Q 2 ) 

2tt 



2tt 

1-Z2 



du 



(14) 



where Q 2 is the momentum scale, a s (Q 2 ) is the strong 
coupling at that momentum scale. 



On the left hand side of Eq. ( 14 ) , the rate that the 
DFFs change with respect to \nQ r ^ is represented. The 
first term on the right hand side represents the effect of 
the parton i emitting a parton j with light-cone momen- 
tum fraction u, with probability Pji(u) that it produces 
the two detected hadrons, hi and while the second 
term represents the effect of two partons, j and k, being 
emitted by i with light-cone momentum fractions u and 
1 — u, respectively, with probability P'j k (u) and each of 
these partons producing one of the detected hadrons. 



r 



We developed a computer code to perform the QCD 
evolution of the dihadron fragmentation functions ac- 



cording to Eq. (14 1, where the DFFs were separated 
into non-singlet, singlet and gluon dihadron fragmenta- 
tion functions. The code is based on the single hadron 
fragmentation function evolution program by the authors 
of Refs. [82 - 84 . The details on the evolution method em- 
ployed, along with the full set of the results will be pre- 
sented in our upcoming paper |85j . As an example, here 
we present the results for LO evolution of D™ + ™ (zi, z%) 
from our model scale of 0.2 GeV 2 to the typical hadronic 
scale of 4 GeV 2 . The results for the evolved DFFs are 
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presented in Figs. 9(a) and |9(b)| corresponding to the 
solutions at z\ = 0.5 and z 2 = 0.5, respectively. The 
dotted red line represents the solution at the model scale 
(0.2 GeV 2 ) and the solid black line represents the solu- 



tion at th e final scale Q 2 = 4 GeV 2 . Both Figs. |9(a) 
and |9(b)| show a shift in the peak of the model results 
towards the lower z region after the evolution, similar to 
the single hadron evolution. 
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(b)z 2 = 0.5: z x D^^ 

FIG. 9. u — ¥ tt + tv~ dihadron fragmentation functions are 
represented as the dotted red lines for the solutions at model 
scale (Qq — 0.2 GeV 2 ) and solid black l ines for the soluti ons 
at the evolved scale (Q 2 = 4 GeV 2 ) for 
z 2 = 0.5. 



VI. CONCLUSIONS 

In this article we have presented results for dihadron 
fragmentation functions calculated within the NJL-jet 
model. DFFs were obtained as numerical solutions of 
the corresponding integral equations derived using the 
quark-jet description of the hadronization process. In 
Section |IV A[ we showed that the integral term, that 
represents the effects of initial undetected hadron emis- 
sion, has a very small effect on the DFFs, except when 



the driving function was zero or when z\ was low. For 
driving functions equal to zero, the corresponding DFFs 
were generated entirely by the integral term and when 
zi was lowered to 0.1, the relative contribution of the 
driving function to the DFF was also lowered for most 
values of z 2 . For the u — > ir + ir~ DFF, the peak value 
of z 2 D^(z u z 2 ) at z x = 0.5 was almost ten times the 
peak value at z\ = 0.1. Because of the lower value of 
the DFF at low z±, the integral term contribution be- 
comes a more significant part of the DFF, reducing the 
relative contribution of the driving function. This effect 
occurs when the fragmenting quark is the favored quark 
for the hadron that receives a small light-cone momentum 
fraction, but is unfavored for the hadron that has access 
to most of the light-cone momentum of the fragment- 
ing quark. One e xamp le where this effect is particularly 
visible is in Fig. 5(c) where the u — > K + K~ DFF is 
mostly composed of the integral term at low z 2 and the 
driving function at higher z 2 . In all three results with 
low z\ (Fig. [5]), the effect is seen for the up quark DFF, 
which is the favored quark for the hadrons h\ — n + and 
hi = K + , but is unfavored for hadrons h 2 — and 
h 2 =K-. 

In Section [lVB| we showed that the strange quarks in- 
clusion has a significant impact on the DFFs. The main 
change to the DFFs when the strange quark is included is 
a considerable reduction in magnitude, similar to the sin- 
gle hadron fragmentations functions [53], caused by the 
availability of the kaon emission channels. The compari- 
son plots for the 7r + 7r~ dihadron fragmentation functions 
were shown in Fig. [6] where the results are similar for the 

7T + 7T + and 7T 



DFFs. 

We examined K 



, — ; , (21,22), where either z\ or z 2 is 

fixed (Figs.[7]and[8|, in Section [TVC] The n+K~ DFFs 
were chosen as they are favored both for a light quark and 
a strange quark. At low values of Z\ and z 2 , the strange 
quark and up quark DFFs, respectively, are dominant. 
The strange quark DFF is dominant for low z\ because 
the 7r + has a small fraction of the light-cone momentum 
of the initial quark, leaving most of the initial momentum 
available to the strange quark's favored fragmentation to 
JV . Similarly, the up quark is dominant for low z 2 be- 
cause the K~ has a small light-cone momentum fraction, 
allowing the the ir + to access most of the momentum for 
its favored fragmentation. Increasing the fixed value of 
either z\ or z 2 to 0.5 in the it + K~ DFF shows that the 
strange quark DFF is dominant. This can be easily inter- 
preted within the model, as only the strange quark can 
produce both the K~ then the ir + in the first two steps 
of the decay chain, whereas the up and down quarks both 
require at least three steps in the decay chain to produce 
both hadrons. 

Finally, in Section [V] we discuss the QCD evolution of 
the dihadron fragmentation functions, which is essential 
in comparing our model calculations with the experimen- 
tal extractions. The evolution equations are presented 
and the method for the numerical solutions is briefly 
discussed. As an example, the results for Df t 77 pre- 
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sented in Fig. [9j show the significant modification of the 
DFFs with evolution. The details for solving the evolu- 
tion equations and the complete set of results for evolved 
DFFs will be presented in our upcoming paper |85j . 
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